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ABSTRACT 

We report the detection of orbital decay in the 12.75-min, detached binary white dwarf (WD) SDSS 
J065133. 338+284423. 37 (hereafter J0651). Our photometric observations over a 13-month baseline 
constrain the orbital period to 765.206543(55) s and indicate the orbit is decreasing at a rate of 
(—9.8 ± 2.8) x 10~ 12 s s _1 (or —0.31 ± 0.09 ms yr _1 ). We revise the system parameters based on 
our new photometric and spectroscopic observations: J0651 contains two WDs with M\ = 0.26 ±0.04 
Mq and M% = 0.50 ± 0.04 M . General relativity predicts orbital decay due to gravitational wave 
radiation of (— 8.2± 1.7) x 10~ 12 s s _1 (or —0.26 ±0.05 ms yr -1 ). Our observed rate of orbital decay 
is consistent with this expectation. J0651 is currently the second- loudest gravitational wave source 
known in the milli-Hcrtz range and the loudest non-interacting binary, which makes it an excellent 
verification source for future missions aimed at directly detecting gravitational waves. Our work 
establishes the feasibility of monitoring this system's orbital period decay at optical wavelengths. 
Subject headings: binaries: close binaries: eclipsing Stars: individual (SDSS 

J065133.338±284423.37) — white dwarfs — gravitational waves 



1. INTRODUCTION 

The 12.75-minute orbital perio d detached b i nary W D 
system J0651 was discovered by iBrown et al.l (|2011l ) as 
part of the extremely low mass (ELM, < 0.25 M ) WD 
Survey, a targeted spectroscopic search for ELM WDs. 
While that survey has yielded some two dozen merger 
systems, with orbit al periods of tens of minutes to hours 
(IBrown et al.llMTot iKilic et al.lf20Tll IBrown et afll20Tl 
iKilic et al.ll2012| ). none are as compact as J0651. 

In addition to a large radial-velocity amplitude, this 
double degenerate system is oriented in such a way that 
it yields a wealth of photometric information: eclipses of 
each star by the other, ellipsoidal variations and Doppler 
boosting. While photometric observations engender an 
accurate way to measure the orbital and system param- 
eters, they also provide multiple clocks with which to 
monitor the orbital evolution of the system. 

The orbital decay of compact binary systems is 
currently the best method to detect the influence 
(and existence) of gravitational waves, and few known 
systems are radiating as strongly or decaying as 
rapidly as J0651. There are presently just five bi- 
naries known with orbital periods less than 15 min- 
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utes, and the other four are interacting: three are 
the AM CVn systems HM Cnc, V407 Vul, and ES 
Cet, and th e other is the l ow-mass X-ray binary 
4U 1820-30 (llsrael et al.l fl999t IHaberl fc Motchl iHfflil 
Steeghs et all 120061: 1 Warner fc WoudtH2002t IStella'etlll 
19871). After~the 5 4- minute HM Cnc (jlsrael et al.1 120021 : 
Roelofs et al.ll2010h . J0651 is the second-loudest gravita- 
tional wave source known in the milli-Hertz frequency 
range (jAmaro-Seoane et al.l I2012T) . J0651 is thus the 
shortest-period detached compact binary known and the 
cleanest system to observe at optical wavelengths for or- 
bital decay due to gravitational wave radiation. 

In this Letter we present follow-up photometric and 
spectroscopic observations of J0651, refine orbital and 
system parameters, and report the detection of rapid or- 
bital decay in the system. We discuss the orbital period 
change in the context of expectations from general rela- 
tivity, as well as deviations expected due to tidal interac- 
tions. Section 2 describes our observations, and Sections 
3 and 4 present the refined system parameters and the 
orbital period decay, respectively. 

2. OBSERVATIONS AND REDUCTIONS 
2.1. Photometric Observations 

Our discovery observations of J065 1 (go = 19.1 mag) 
were described bv lBrown et al.l (j 2 1 lh and included some 
12.7 hr of photometry from the McDonald Observatory 
(McD) 2.1 m Otto Str uve Telescope using the A rgos 
frame-transfer camera dNather fc Mukadam 2004]). In 
the subsequent year, we have obtained an additional 
195.4 hr of photometry using four different facilities. 

The majority of our continuing photometric observa- 
tions were carried out in an identical manner to the dis- 
covery observations, with the Argos instrument through 
a 1 mm BG40 filter to reduce sky noise. We obtained 
5 — 30 s exposures of the target, with a typical exposure 
time of 10 s, depending on the observing conditions. Ad- 
ditionally, we obtained nearly 3 hr of data in December 
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Fig. 1. — Spectroscopic observations of J0651. The top left panel 
shows the summed spectra from our October 2011 and April 2012 
observations, with model fits to the H/3 to Hll lines, which we 
use to derive the primary parameters in Section [3] The top right 
panel shows fits to the lower S/N, summed spectra from March 
2011 from [Brown ct al. (201lj). The middle panel shows our new 
radial velocity observations of J0651 from three new epochs, and 
the bottom panel shows those data phased to the orbital period. 



2011 using the Agile instrument (jMukadam et al.ll2011[) 
mounted on the 3.5 m telescope at Apache Point Obser- 
vatory (APO), using a 1 mm BG40 filter and 10 — 15 s 
exposures. In January and March 2012, we obtained 6.8 
hr of data with 20 s exposures using GMOS (IHook et al.1 
I2004T ) on the 8.1 m Gemini North telescope as part of 
the queue programs GN-2011B-Q-95 and GN-2012A-Q- 
29. Most of the Gemini data were taken using a Sloan-g 
filter, but we obtained nearly 2 hr of data using a Sloan-r 
filter to constrain the luminosity of the secondary WD 
(see Section [3} . Additionally, we obtained 1.5 hr of data 
using 10 s exposures in March 2012 and 1.0 hr of data 
using 5 s exposures in Apr i l 2012 using the OSIRIS in- 
strument (|Cepa et al.ll2000l [200l through a Sloan-g fil- 
ter and in fast photometry mode, mounted on the Gran 
Telescopio Canarias (GTC) 10.4 m telescope. 

We bias- and flat-field correct the raw science frames 
using standard IRAF routines. For Argos and Agile, 
we perform weighted, circular, fixed-aperture photom- 
etry on the calibrated frames using the external IRAF 



package ccdJisp ()Kanaan et al.l I2002D . We divide the 
sky-subtracted light curves using five brighter compar- 
ison stars in the field to remove transparency varia- 
tions. To remove any long-term trends caused by differ- 
ential atmospheric extinction, we fit a low-order polyno- 
mial to observi ng runs exceeding 2 hr using the WQED 
software suite ([Thompson fc Mullally 2009), which we 
also use to apply a timing correction to each obser- 
vation to account for the motion of th e Earth aroun d 
the barycenter of the solar system (iStumpfB Il980f). 
We u se the formalism described in lEverett fe Howell 
(|2001D to derive average point-by-point photometric er- 
rors of 1.0 mmag for GMOS and OSIRIS, and 2.8 mmag 
for Argos observations. We calibrate these errors us- 
ing the go = 19.1 mag, photometrically constant star 
SDSS J065132. 86+284408.4, within 20" of our target. 

For the GM OS and OSIRIS data, we use DAOPHOT 
([Stetson! [1987) to perform aperture photometry on our 
target and a dozen photometrically constant SDSS point 
sources in our images for calibration. We use a script by 
lEastman et al.l ([2010D to apply a barycentric timing cor- 
rection and check it against the timings based on WQED. 

2.2. Spectroscopic Observations 

To obtain higher signal-to-noise-ratio spectroscopy and 
better phase coverage, we obtained additional time-series 
spectroscopy of J0651 at the 6.5m MMT telescope in Oc- 
tober 2011 and April 2012. All 79 sp ectra were taken us- 
ing the Blue Channel Spectrograph ([Schmidt et al.lfl989h 
with a 800 1 mm -1 grating and a 1" slit. This set-up 
provides 2.1 A resolution and a spectral coverage from 
3550 — 5450 A. The lower resoluti on compared to the dis- 
covery observations presented bv I Brown et al.l ([201 1! ) en- 
abled us to decrease the exposure time to 90 s, providing 
a radial velocity curve with better phase coverage (see 
Figure [J). The reduced spectra have an average radial 
velocity error of ±31 km s _1 . 

3. UPDATED SYSTEM PARAMETERS 

We compute th e orbi tal elements using the code of 
iKenyon fc Garcia! (|1986| ). which weights each velocity 
measurements by its associated error. However, the ob- 
served velocity amplitude is an underestimate because 
our 90 s exposures span 12% of the orbital phase, and 
because the radial velocity curve is not linear. By inte- 
grating a sine curve at the phase covered by our expo- 
sures, we determine that the velocity amplitude correc- 
tion is 2.3%. The resultant best-fit, corrected velocity 
semi-amplitude is K = 616.9 ± 3.5 km s _1 . This is sig- 
nificantly lower than our original value, 657.3 ± 2.4 km 
s _1 , computed in the same way. 

Given that the original radial velocity semi-amplitude 
depended on a dozen measurements at quadrature with 
30 km s _1 errors, the formal error was an underestimate. 
A Monte Carlo calculation indicates that the true uncer- 
tainty in K using our original 27 spectra was actually 14 
km s _1 . The same calculation with the new data yields 
an uncertainty in if of 5 km s _1 . This implies that the 
companion i s less massive than originally predicted by 
iBrownet all (p?5ll . 

We refine the physical parameters of the primarjQ WD 



1 Following IjBrown et alJ[20TIh we refer to the low-mass WD as 
the primary since it contributes >95% of optical light. 
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Fig. 2. — The top panels show high-speed photometry of J0651 from the 8.1 m Gemini North telescope and the GTC 10.4 m telescope, 
folded at the orbital period and duplicated for clarity. Directly below is the same data binned into 100 phase bins, with error bars, and 
over-plotted with our best-fit model. The bottom panel shows a portion of the folded, binned data from four different subsets, a visual 
representation of the (O — C) diagram in Figure [3] The decreasing orbital period is evident as the primary eclipses, shown, shift sooner. 



using the summed spectrum, which has S/N=78 per res- 
olution element. Fitting our new spectra with the ste l- 
lar atmosphere models of iTremblav fc Bergeron! ([2009) , 
which include improved Stark broadening profiles with 
non-ideal gas effects, formally yields T a g = 16530 ± f 05 
K and logg = 6.76 ± 0.02. This result is nearly identical 
to our original mea surements (Tiff = 1 6400 ± 300 K and 
logg = 6.79 ±0.04. [Brown et alBMTl ). 

Additionally, we investigate the effect of velocity 
smearing on the derived atmospheric parameters by ana- 
lyzing the spectra obtained at quadrature (\v\ > 500 km 
s _1 , when velocity smearing should be at its minimum). 
For these spectra at quadrature, we find that T e g is 500 
K lower and logg is 0.07 dex higher. These differences 
reflect our systematic error, and also indicate how the 
parameters of the tidally distorted primary depend on 
phase, but the parameters remain consistent with the 
higher S/N summed and phased spectrum. Thus we 
adopt a mean T cff = 16530±20 K and log q = 6. 76±0.04, 
implying a 0.25 M Q primary (|Panei et al.ll2007l ). 

We use our high-quality g-band Gemini and GTC pho- 



tometry to refine the system p arameters using the ligh t 
curve fitting code JKTEBOP (|Southworth et al.|[200l . 
We supply the limb-darkening coefficients from WD at- 
mosphere models appropriate for the J0651 system using 
I ((/,)/ 1(1) = 1 — ci(l — fx) — €2(1 — y/JJ), where fi — cos6> 
(P. Bergeron 2012, private communication). These co- 
efficients are included in Table [TJ their uncertainties are 
negligible given our observed T e g and logg uncertain- 
ties. Additionally, we have adopted gravity-darkening 
coefficients of (3i = /?2 — 0.36 for both the primary and 
secondary, where F oc T^ ff oc g^. We expect convection 
to be present in both stars, and our light curve fits do 
poorly for ft = 1.0, as expected for a purely radiative 
atmosphere, so adopting /3 = 0.36 is reasonable. 

We fix the limb- and gravity-darkening coefficients and 
fit for the inclination and component radii, and our error 
estimates resul t from 10,000 Mo n te Ca rlo simulations, 
as described in ISouthworth et al.l (|2005l ). Gravitational 
lcnsing should minimally affect the derived radius of the 
primary a nd secondary by roughly 0.1% and 0.7%, re- 
spectively (Marsh 2001), and have not been included in 



4 



Hermes et al. 



TABLE 1 




System parameters. 




Parameter 


Value 


[method used to derive parameter] 





TABLE 2 
Journal of mid-eclipse times. 



Orbital Period [phot.] 

K\ (corrected for smearing) [spec] 

7„ ei [spec] 

Primary T c ff [spec] 

Primary log g [spec] 

Primary Mass (Mi) [phot.] 

Primary Radius (-Ri) [phot.] 

Inclination (i) [phot.] 

Mass Ratio (q) [phot., spec] 

Secondary Mass (M2) [spec] 

Secondary T G g [phot.] 

Secondary Radius (R2) [phot.] 

Limb Darkening, Primary, g-band 

Limb Darkening, Secondary, g-band 

Limb Darkening, Primary, r-band 

Limb Darkening, Secondary, r-band 



765.206543(55) s 
616.9 ±5.0 km s" 1 
-7.7 ±4.5 km s" 1 
16530 ± 200 K 
6.76 ±0.04 
0.26 ± 0.04 M 
0.0371 ± 0.0012 R @ 
84.4 ± 2.3 degrees 
1.92 ±0.46 
0.50 ± 0.04 M© 
8700 ± 500 K 
0.0142 ± 0.0010 R e 
Cl = -0.106, c 2 = 0.730 
ci = -0.128, c 2 = 0.898 
Cl = -0.076, c 2 = 0.562 
ci = -0.099, c 2 = 0.735 



the fits. The primary radius is a volume-average; the 
tidal distortions make the star 3.3% oblate. 

The photometry allows us to test the ELM WD mod- 
els by providing an independent estimate on the mass of 
the primary. To do so, we hold fixed a series of differ- 
ent mass ratios in our light curve fits, and in each case 
use the resultant secondary radiu s in combina tion with 
the tested mass-radius relation of iWoodl (|1995f ) in order 
to back out the mass of the primary. Consistently, this 
method finds M x = 0.26 ± 0.04 M Q , which we adopt. 

Taking Ah = 0.26 ± 0.04 M© and K = 616.9 ± 5.0 km 
s , the secondary mass is thus M~ = 0.50±0.04 M© for 
the best-fit inclination of 84.4 ± 2.3 deg. Table Q] shows 
our final light curve results found by fixing q — 1.92. 
The resulting radius of the secondary P 2 = 0-0142 ± 
0.0010 i? Q , implies M 2 = 0.50 ± 0.04 M© (IWoodl 119951 ). 
in good agreement. Pairing the volume-averaged primary 
radius R\ = 0.0371±0.0012 P© with the observed surface 
gravity yields Mi = 0.29 ± 0.05 M©, somewhat larger 
but consistent with our adopte d valu e, as well as with 
the result using the iPanei et al.l (|2007|) models. 

Finally, we use our Gemini r-band data to constrain 
the luminosity and temperature of the secondary. Fixing 
the limb- and gravity-darkening coefficients, and adopt- 
ing the inclination and component radii from the g-band 
fits, we find the secondary contributes 3.7 ±0.2% of light 
in the g-band and 4.6 ± 0.6% of light in the r-band. 

Adopting M g = 8.9 ±0.1 ma g and M r = 9.2 ± 0.1 
mag for the 0.26 M© primary (jPanei et al.l I2007T ) . the 
secondary thus has M g = M r = 12.5 ± 0.2 mag. For a 
0.5 M© WD, cooling modek0 suggest a temperature of 
8700 ± 500 K for the secondary, which corresponds to 
a coo l ing age of roughly 700 My r dHolberg fc Bergeron! 
20061: iKowalski fc Saumonl 120061 : iTremblav et al.l (20111 : 
Berg eron et al.H201lD .~ 



4. DETECTION OF ORBITAL PERIOD DECAY 

We demonstrate a secular change in the orbital period 
of J0651 by constructing an (O — C) diagram, where we 
compare the observed mid-eclipse times (O) to expected 
mid-eclipse times computed from the assumption of a 
fixed orbital period (C) for future epochs (E — t/P). To 

3 http://www.astro.umontreal.ca/~bergeron/CoolingModels 



Facility 


Epoch 


No. Eclipses 


Mid-Eclipse Time 




(E) 




(BJD TDB - 2450000) 


McD 


373 


59 


5655.9015854(85) 


McD 


20735 


11 


5836.2388047(83) 


McD 


23327 


52 


5859.1949883(81) 


McD 


27251 


146 


5893.9480995(57) 


APO,McD 


30981 


110 


5926.9830564(67) 


Gem-N 


31207 


5 


5928.984606 (11) 


McD 


34164 


284 


5955.1734668(37) 


McD 


36291 


80 


5974.0113686(59) 


GTC 


39171 


7 


5999.518212 (12) 


McD 


39383 


53 


6001.3958308(75) 


Gem-N 


39426 


6 


6001.776617 (14) 


Gem-N 


39542 


10 


6002.803962 (10) 


McD 


43317 


89 


6036.2375115(72) 


GTC 


43446 


5 


6037.3800048(72) 


McD 


46578 


16 


6065.118768 (26) 



Note. — While all four instruments we have used are conditioned 
with a GPS receiver and thus should have absolute time stamps ac- 
curate to a few ms, the mid-exposure times for our GMOS-N data 
all end in .2 or .7 s, suggesting a rounding error of up to 0.25 s; the 
uncertainties in our Gemini mid-cclipsc times have thus been enlarged 
by 0.25 s. We may remove any potential systematic time offsets by 
computing the period change using only the points from McDonald 
Observatory, which yields (-8.2 ± 3.2) X 10 -12 s s" 1 . 

estimate the mid-eclipse times, we fix the best-fit model 
parameters from our analysis in Section [3] and fit each 
subset of observations only for the mid-eclipse time near- 
est the mea n time of the observ ations. 

Following Kepler et al.l (|1991l ). if the orbital period is 
changing slowly with time, we can expand the observed 
mid-time of the E th eclipse, tg, in a Taylor series around 
Eg to arrive at the classic (O — C) equation 



O - C = AT, 



AP E + ^P PE 2 



where T is the mid-time of the first eclipse, AT- is the 
uncertainty in this mid-point, Po is the orbital period at 
the first eclipse and AP is the uncertainty in this period. 
Any secular change in the period, dP/dt, will cause a 
parabolic curvature in an [O — C) diagram. Currently, 
the acceleration in the period change, d(dP/dt)/dt, is 
negligible, and we will limit our discussion to a second- 
order polynomial fit. 

To construct an (O — C) diagram, we must first de- 
termine To and Po- A preliminary estimate comes from 
a simple Fourier transform of our whole data set, which 
we use to create an initial (O — C) diagram. We then it- 
eratively adjust T and Po by the zeroth- and first-order 
terms from our best-fit parabola until the adjustments 
are smaller than the error in these terms; these errors re- 
sult from the covariance matrix. Our recomputed, final 
(O — C) diagram uses this new ephemeris and period and 
is shown in Figure [3l We find: 

T = 2455652.5980910 ± 0.0000084 BJD T db 

P = 765.206543 ± 0.000055 s 

Table [2] presents the mid-eclipse times from each subset 
of our observations. Each night of observing from Gemini 
and GTC have been given their own subset, as have each 
month of data from McDonald. Since AXb and APo are 
nonzero, the zeroth- and first-order terms of the parabola 
indicating the predicted dP/ dt in Figure [3] have been 
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Fig. 3. — (O — C) diagrams of the orbital evolution in J0651 since April 2011; blue dots represent data from McDonald Observatory 
and APO, green squares from Gemini-North, and maroon triangles from GTC. The top panel shows the change in mid-eclipse times as 
determined by light curve modeling, and the best-fit parabola yields an estimate for the observed rate of orbital period change. Additionally, 
the bottom panel shows the results from a model-independent, linear least-squares fit using the orbital period and higher harmonics. The 
dotted line at (O — C) = shows the line of zero orbital decay, while the grey dashed line shows the predicted orbital decay expected solely 
from gravitational wave radiation. Using both methods, our early results match the GR prediction to the \-a level. 



allowed to vary within the current constraints on these 
terms. 

A weighted, second-order, least-squares fit to the mid- 
eclipse times yields a rate of period change of (—9.8 ± 
2.8) x 1CT 12 s s- 1 (or -0.31 ±0.09 ms yr" 1 ). This value 
includes our May 20f 2 data point, which has just 3.4 hr 
of data spread over four nights at a minimum airmass of 
2.0. If we do not include this last point, the inferred rate 
of period change differs slightly, yielding (— 10.6 ± 2.9) x 
10~ 12 s s -1 (or -0.33 ± 0.09 ms yr" 1 ). A parabola is 
needed to best represent the data: The best second-order 
fit has x 2 = 33.0 (12 d.o.f.), whereas the best first-order 
fit has x 2 = 44.9 (13 d.o.f.). 

As a sanity check, we also construct an (O—C) diagram 
using a model-independent approach. Here we perform a 
simultaneous least-squares fit to each subset of data using 
a series of sine curves at the orbital period (Po) up to the 
last harmonic before the Nyquist frequency of that sub- 
set. This effectively uses the high-amplitude ellipsoidal 
variations at half the orbital period as our clock, with a 
Fourier series of the orbital harmonics to reproduce the 
eclipses. A best-fit parabola to the observed minima of 
the ellipsoidal variations is shown in the bottom panel of 
Figure[3]and yields dP/dt = (-5.1 ± 1.8) x 10~ 12 s s _1 . 
In order to compare this to the mid-eclipse times at the 



orbital period, we must multiply this result by a factor 
of two, which yields dP/dt = (-10.2 ±3.5) x 10~ 12 s s" 1 
(or -0.32 ± 0.11 ms yr" 1 ). 

This model- independent (O — C) method retains larger 
errors because the harmonics are inevitably truncated 
by the Nyquist frequency in the observations and so 
are only roughly capable of replicating the deep primary 
eclipses. Our results from this method are not orthogo- 
nal to the model-dependent light curve fitting, as both 
fit the eclipses and ellipsoidal variations. Therefore, both 
results cannot be averaged, and we emphasize the results 
from the mo del- dependent approach. 

Thus, our best estimate for the rate of orbital period 
change in J0651 after 13 months is (-9.8 ± 2.8) x 10" 12 
s s _1 (or —0.31 ± 0.09 ms yr -1 ), a 3-cr detection. This 

yields a timescale for period change, P/P — 2.5 ± 0.8 
Myr. 

5. DISCUSSION 

Based on the refined parameters for the J0651 system 
(see Section[3]) and treating each WD as a point mass in a 
non-relativistic circular orbit, general relativity predicts 
an orbital period decay in this system of (—8.2 ± 1.7) x 
10~ 1 2 s s -1 (or -0. 26±0.05ms vr" 1 ) (lLandau fc Lifshitzl 
Il975f) . Recently, Ivan den Broek et alj (|2012D demon- 
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strated that the point mass approximation is valid to 
better than 1% in cases such as J0651, so our uncer- 
tainty in the orbital decay from gravitational wave ra- 
diation is dominated by the uncertainty in the compo- 
nent masses. Additional effects that could modulate the 
mid-eclipse times are unlikely to explain the observed 
shift: For example, give n the estimated distance of 1.0 
kpc (jBrown et al.ll201ll) . we expect the proper motion 
to change the pe riod by no more than 5 x 10 -16 s s _1 
(IShklovskiilll970h . 

It is evident from the high-amplitude ellipsoidal vari- 
ations of the primary that strong tidal forces are also 
present. These tides will act as a torque to spin- up the 
WDs if the system is synchronized, further robbing the 
orbit of angular momentum and increasing the rate of 
orbital period decay. The degree to which this tidal 
torquing influences the orbital evolution depends on the 
effective tidal locking, which is in many ways determined 
by the physical structure of the ELM WD. This effect 
could increase the rate of p eriod decay by at least 5% if 
the system is synch ronized (jPirdl20 1 Tl lBenacquistall201 li 
iFuller feT ai 2012). With just 13 months of monitoring, 
our sensitivity in the observed rate of orbital decay is 
not yet sufficient to detect a significant deviation from 
pure gravitational wave losses. However, future obser- 
vations should constrain this discrepancy, providing an 
excellent probe of the interior of ELM WDs, in addi- 
tion to the possibilities opened through asteroseismology 



(jSteinfadt et al.ll20lol : iHermes et al.|[2012l ). 

The short period of J0651 makes it one of the loudest 
known sources of gravitational wave radiation, and con- 
tinued monitoring of orbital decay in the system will pro- 
vide strong constraints on the gravitational wave strain 
of J0651. This is important for future gravitational wave 
missions, like the evolved Laser Interferometry Space An- 
tenna (eLISA). Critical to that success is disentangling 
the contributions of tidal torques on the orbital decay, 
an effort worthy of further photometric observations and 
modeling. 
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rold. C.A.P. is thankful to Alessandro Ederoclite for 
suggesting to observe J0651 with GTC, and to Sebastian 
Hidalgo for help with photometric reductions. J.J.H., 
M.H.M. and D.E.W. acknowledge the support of the NSF 
under grant AST-0909107 and the Norman Hackerman 
Advanced Research Program under grant 003658-0252- 
2009. A.S.M. acknowledges observing support from NSF 
grant AST-1008734. 

Facilities: McDonald 2.1m (Argos), MMT (Blue 
Channel Spectrograph), Gemini North (GMOS), GTC 
(OSIRIS), APO (Agile) 



REFERENCES 



Amaro-Seoane, P., Aoudia, S., Babak, S., et al. 2012, 

larXiv:1201.3621l 
Bcnacquista, M. J. 2011, ApJ, 740, L54 

Bergeron, P., Wesemael, F., Dufour, P., et al. 2011, ApJ, 737, 28 
Brown, W. R., Kilic, M., Allende Prieto, C, & Kenyon, S. J. 

2010, ApJ, 723, 1072 
Brown, W. R., Kilic, M., Hermes, J. J., et al. 2011, ApJ, 737, L23 
Brown, W. R., Kilic, M., Allende Prieto, C, & Kenyon, S. J. 

2012, ApJ, 744, 142 
Cepa, J., Aguiar, M., Escalera, V. G., et al. 2000, Proc. SPIE, 

4008, 623 

Cepa, J., Aguiar-Gonzalez, M., Bland-Hawthorn, J., et al. 2003, 

Proc. SPIE, 4841, 1739 
Eastman, J., Siverd, R., & Gaudi, B. S. 2010, PASP, 122, 935 
Everett, M. E., & Howell, S. B. 2001, PASP, 113, 1428 
Fuller, J., & Lai, D. 2012, MNRAS, 421, 426 
Haberl, F., & Motch, C. 1995, A&A, 297, L37 
Hermes, J. J., Montgomery, M. H., Winget, D. E., et al. 2012, 

ApJ, 750, L28 
Holberg, J. B., & Bergeron, P. 2006, AJ, 132, 1221 
Hook, I. M., j0rgensen, I., Allington-Smith, J. R., et al. 2004, 

PASP, 116, 425 

Israel, G. L., Panzera, M. R., Campana, S., et al. 1999, A&A, 
349, LI 

Israel, G. L., Hummel, W., Covino, S., et al. 2002, A&A, 386, L13 
Kanaan, A., Kepler, S. O., & Winget, D. E. 2002, A&A, 389, 896 
Kenyon, S. J., & Garcia, M. R. 1986, AJ, 91, 125 
Kepler, S. O., Winget, D. E., Nather, R. E., et al. 1991, ApJ, 378, 
L45 

Kilic, M., Brown, W. R., Allende Prieto, C, et al. 2011, ApJ, 
727, 3 

Kilic, M., Brown, W. R., Allende Prieto, C, et al. 2012, ApJ, 

751, 141 

Koester, D. 2008, arXiv:0812.0482 

Kowalski, P. M., & Saumon, D. 2006, ApJ, 651, L137 

Landau, L. D., & Lifshitz, E. M. 1975, The Classical Theory of 

Fields. Pergamon Press, Oxford 



Marsh, T. R. 2001, MNRAS, 324, 547 

Mukadam, A. S., Owen, R., Mannery, E., et al. 2011, PASP, 123, 
1423 

Nather, R. E. & Mukadam, A. S. 2004, ApJ, 605, 846 

Panei, J. A., Althaus, L. G., Chen, X., & Han, Z. 2007, MNRAS, 

382, 779 
Piro, A. L. 2011, ApJ, 740, L53 

Roelofs, G. H. A., Rau, A., Marsh, T. R., et al. 2010, ApJ, 711, 
L138 

Schmidt, G. D., Weymann, R. J., & Foltz, C. B. 1989, PASP, 
101, 713 

Shklovskii, I. S. 1970, Soviet Ast., 13, 562 

Southworth, J., Maxted, P. F. L., & Smalley, B. 2004, MNRAS, 
351, 1277 

Southworth, J., Smalley, B., Maxted, P. F. L., Claret, A., & 

Etzel, P. B. 2005, MNRAS, 363, 529 
Steeghs, D., Marsh, T. R., Barros, S. C. C, et al. 2006, ApJ, 649, 

382 

Steinfadt, J. D. R., Bildsten, L., & Arras, P. 2010, ApJ, 718, 441 
Stella, L., Priedhorsky, W., & White, N. E. 1987, ApJ, 312, L17 
Stetson, P. B. 1987, PASP, 99, 191 
Stumpff, P. 1980, A&AS, 41, 1 

Thompson, S. E., & Mullally, F. 2009, Journal of Physics 

Conference Series, 172, 012081 
Tremblay, P.-E., & Bergeron, P. 2009, ApJ, 696, 1755 
Tremblay, P.-E., Bergeron, P., & Gianninas, A. 2011, ApJ, 730, 

128 

van den Broek, D., Nclcmans, G., Dan, M., & Rosswog, S. 2012, 

larXiv: 1206.07441 
Warner, B., & Woudt, P. A. 2002, PASP, 114, 129 
Wood, M. A. 1995, in Proc. 9th European Workshop on White 

Dwarfs, ed. D. Koester & K. Werner (Berlin: Springer), 41 



